Ground states of a frustrated spin-^ antiferromagnet: Cs 2 CuCl 4 in a magnetic field 
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We present detailed calculations of the magnetic ground state properties of CS2CUCI4 in an applied magnetic 
field, and compare our results with recent experiments. The material is described by a spin Hamiltonian, deter- 
mined with precision in high field measurements, in which the main interaction is antiferromagnetic Heisenberg 
exchange between neighboring spins on an anisotropic triangular lattice. An additional, weak Dzyaloshinkii- 
Moriya interaction introduces easy-plane anisotropy, so that behavior is different for transverse and longitudinal 
field directions. We determine the phase diagram as a function of field strength for both field directions at zero 
temperature, using a classical approximation as a first step. Building on this, we calculate the effect of quantum 
fluctuations on the ordering wavevector and components of the ordered moments, using both linear spinwave 
theory and a mapping to a Bose gas which gives exact results when the magnetization is almost saturated. Many 
aspects of the experimental data are well accounted for by this approach. 

PACS numbers: 75.10.Jm, 75.25.+z,75.40.Gb 



I. INTRODUCTION 

The layered, insulating magnet CS2C11CI4 has at- 
tracted intense recent experimental and theoretical 
attentioni 1 i 2 i 3 i 4 i 5 i 6 i 7 i 8 i 9 i '"i 1 1 Much of the interest arises be- 
cause the material is a rare example of a spin S = |, 
triangular lattice antiferromagneti^i Its small spin, quasi- 
two dimensionality and geometric frustration are all features 
expected to enhance zero-point fluctuations in Neel-ordered 
states, and to promote spin-liquid states»i2iiii4 Indeed, 
inelastic neutron scattering experiments on CS2CUCI4 4 
have revealed extended scattering continua in the dynamic 
structure factor, and various spin liquid states 7 ' 911 have been 
proposed to explain this observation in terms of fractionalized 
excitations. Nevertheless, at sufficiently low temperature 
CS2CUCI4 displays conventional, magnetically ordered states 
over much of the phase diagram spanned by magnetic 
field strength and directioniiSii 4 . In this paper, we develop 
a theoretical treatment of these ordered states using two 
approaches. Starting from a classical description, we discuss 
fluctuations using linear spinwave theory. And starting from 
the fully polarized spin state reached at high field, we discuss 
fluctuations as a dilute Bose gas of spin flips. We present a 
detailed comparison of our results with experiment. 

The low-temperature states of CS2CUCI4 have been exam- 
ined as a function of magnetic field strength and direction, 
using neutron diffraction. 2 ' 3 At zero field, long-range order, in 
the form of an incommensurate spiral spin structure, occurs 
below a Neel temperature of T/v = 0.62 K. The magnetic mo- 
ments lie in an easy plane due to anisotropy arising from a 
Dzyaloshinskii-Moriya interaction. The presence of this in- 
teraction, breaking SU(2) symmetry in the spin Hamiltonian, 
has profound consequences for the behavior of the system in a 
magnetic field, and the ordering observed depends on the field 
direction. 2 Two field directions have been studied: transverse 



to the easy plane (along the crystallographic a direction) and 
longitudinal - within the easy plane (the crystallographic b-c 
plane). In a transverse field, spins cant out of the easy plane 
towards the field direction, gaining Zeeman energy. Below a 
critical field of B" r = 8.44 T, ordered moments at different 
sites lie on a cone around the field direction. As the critical 
field is approached the cone angle closes to zero, and above it 
the magnetization is saturated. The behavior for a longitudi- 
nal magnetic field is considerably more complex^ For fields 
along the c axis of strength B c , at weak fields, B c < 1.4 T, 
the anisotropy confines the spins in the b-c plane, creating a 
distorted cycloid. In the field range 1.4T < B c < 2. IT, a 
second incommensurately ordered phase appears. At interme- 
diate field strengths, in the range 2.1 T < B c < 7.1 T, no 
magnetic Bragg peaks have so far been reported. In stronger 
fields, magnetic Bragg peaks 5 at incommensurate wavevector 
are found up to the critical field B c cr . = 8.0 T, beyond which 
the magnetization is saturated at low temperature. 

In this article our starting point is the spin Hamiltonian 
for CS2CUCI4, as determined by high-field experiments^ We 
discuss the symmetry of this Hamiltonian and establish its 
ground-state phase diagram in transverse and longitudinal 
fields, within a classical approximation. We find incommen- 
surate phases of three types. Extending our treatment to in- 
clude quantum fluctuations, we proceed in two ways. First, 
we set out linear spinwave theory, treating fluctuations around 
the classical state at leading order using a standard 1/S ex- 
pansion. Although the expansion parameter is not small 
in the case at hand, results known for the nearest neighbor 
Heisenberg antiferromagnet on the squar e) 15 i 16 and isotropic 
triangular 17 lattices suggest that linear spinwave theory is 
likely to be quite accurate even for S = h. Second, supple- 
menting the 1 / S expansion, we apply theory for a dilute Bose 
gas to spin flips in a system with almost saturated magnetiza- 
tion. Using both methods, we determine quantum corrections 
to the ordering wavevector and components of the local or- 



dered moments as a function of field strength. The results de- 
pend markedly on the presence of a Dzyaloshinskii-Moriya 
interaction and on the orientation of the applied magnetic 
field. We also investigate the effect of interlayer exchange, 
focussing on its influence on magnetic order in a transverse 
field. We compare our results extensively with experimental 
data. 

The remainder of the paper is organized as follows. We 
introduce the spin Hamiltonian and discuss its symmetries in 
Sec. Ql] In Sec. [TO] we develop classical theory and estab- 
lish the phase diagram within a classical approximation, in 
transverse and longitudinal fields. We examine the effects of 
quantum fluctuations using the 1/5 expansion and dilute Bose 
gas methods in Sec. IIVI calculating static properties and com- 
paring these with experimental data. In Sec. [V]\ve consider 
interlayer coupling. Finally, in Sec. IV II we summarize our 
conclusions. 

II. CRYSTAL STRUCTURE AND SPIN HAMILTONIAN 
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FIG. 1 : Magnetic sites and exchange couplings within a single layer 
of CS2C11CI4. Layers are stacked along the a-direction, with inter- 
layer spacing a/2 and a relative displacement in the c-direction. 



The magnetic moments in CS2CUCI4 are carried by Cu 2+ 
ions. The orthorhombic unit cell contains four CuCl 2- tetra- 
hedra arranged in two layers in the b-c plane. 4 The location 
of magnetic sites within a single layer is illustrated in Fig. [2 
Exchange interactions are sufficiently weak that it is possible 
using laboratory magnetic fields to fully polarize the moments 
at low temperature, and the spin Hamiltonian has been deter- 
mined from a study of the excitation spectrum in a saturat- 
ing transverse fields This method has the advantage of yield- 
ing interaction constants with the minimum of theoretical as- 
sumptions, since it focusses on the dynamics of single spin 
flips. In this way it has been established that the largest inter- 
action is antiferromagnetic exchange J, coupling neighboring 
spins along the chains, and that neighbors on adjacent chains 
have a weaker exchange coupling J'. In addition, the mea- 
surements indicate a Dzyaloshinskii-Moriya (DM) exchange 
between chains, allowed by symmetry, 18 and a weak antifer- 
romagnetic nearest-neighbor interlayer coupling J", which 
stabilizes long-range magnetic order against thermal fluctu- 
ations. 

The model Hamiltonian, with experimentally determined 
parameter values given in Table|IJ is 

W = Tia + Ti-DM + 7~Lb j (1) 

where Ho is the Heisenberg exchange energy, TLdm repre- 
sents the DM interaction, and Hb is the Zeeman energy in 
an applied magnetic field. Denoting spin-i operators at the 
sites R of a stacked anisotropic triangular lattice by Sr, the 
exchange energy is 

7~i<3 = ^ [-/Sr. • Sr + 6 1 _|_6 2 
R 

+J (Sr • Sr + 5 x + Sr • Sr +< 5 2 ) 

+ J"Sr • Sr +< 5 3 ] , (2) 



TABLE I: Hamiltonian parameters, from Ref. 3. 



Parameters 


Exp. 


J (meV) 


0.374(5) 


J' (meV) 


0.128(5) 


J" (meV) 


0.017(2) 


D (meV) 


0.020(2) 



where the nearest neighbor vectors Si and £2 are indicated 
in Fig. ^ an d the out-of-plane vector ($3 connects spins on 
adjacent layers. The DM energy is 

H DM = - $](-l) n D • Sr x [S R+5l + S R+<52 ] , (3) 

R 

where D = (D, 0, 0) is a vector associated with the oriented 
bond between the two coupled spins and n is a layer index. 
The factor (—1)™ indicates that the interaction alternates be- 
tween even and odd layers, which are inverted version of one 
another. The Zeeman energy arising from a magnetic field 
B = (B a ,B b , B c ) is 

Hb = -Y.9^bB 1 S 1 r , (4) 

R 

where g is the gyromagnetic tensor g = (2.20, 2.08, 2.30). 19 
We omit the dipole-dipole interaction and several small ef- 
fects, including a relative offset of the Cu ions along c between 
adjacent layers, a small component of the D vector perpen- 
dicular to the a-axis and possible anisotropy of the exchange 
interactions in spin space. 

At the classical level, the intrachain coupling J favors a 
staggered magnetization in the spin chains and the interchain 
coupling J' frustrates this state. As J' /J is varied, 7io in- 
terpolates between the fully frustrated Hamiltonian for the 
isotropic triangular lattice ( J' = J) and that for uncoupled 
one dimensional spin chains (J 7 = 0). The DM interaction 
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favors states in which spins lie in the b-c plane, with a rotation 
of 7r/2 between adjacent spin chains. 

It is convenient at this point to introduce notation associated 
with reciprocal space. We express wavevectors in terms of the 
reciprocal lattice vectors, writing Q = (h, k, I) as shorthand 
for 2n(h/a : k/b,l/c). The Fourier transforms of the exchange 
and DM interactions are 



and 



Jq = J cos(2-7rfc) + 2 J' cos(7rfc) cos(7r/) (5) 



D Q = -2D sin(Trfc) cos(ttZ) . (6) 



When considering transverse magnetic fields, these appear in 
the combination 



Jq = Jq+ Dq 



(7) 



We close this section with a discussion of the symmetry 
of the Hamiltonian H. While Ho has full SU(2) spin sym- 
metry, Hdm has a lower, Z2 <g> U(l) symmetry. Here, U(l) 
arises from spin rotations around the D vector, and Z2 orig- 
inates from invariance under the combination of space inver- 
sion (R — > — R) and the spin operations 



Sxi 
S-x 



-S x x 
S-x 



(8) 



where x is an arbitrary unit vector in the b-c plane. To il- 
lustrate the nature of the Z2 symmetry, one can consider the 
chiral scalar K = J2a ' (S2 x S3), where the spin product 
is performed in a cyclical fashion over all triangular plaque- 
ttes. Under the Z2 operation, Ho + Hdm is invariant but 
K — > —K. The inclusion of Hb further reduces the symme- 
try, to U( 1 ) in a transverse magnetic field (with S a a conserved 
quantity), and to Z 2 in a longitudinal field. 



in. CLASSICAL ANALYSIS 

The classical approximation consists of treating the spin op- 
erators S as classical vectors of length S = 1/2. The Hamil- 
tonian then becomes an energy functional which can be min- 
imized to determine the magnetic structure. Omitting inter- 
layer exchange and DM interactions, the classical ground state 
in zero field is a spin spiral 





S R = S \ cos(Q* ■ R + a) 
,sin(Q*; -R + a) 



(9) 



where the arbitrary phase a reflects spontaneous breaking of 
the U(l) symmetry and the wavevector Q* ; is determined 
by minimizing the exchange energy Jq. We find Q* ; = 
±(0,l/2 + e* ,0) where e* = vr" 1 arcsin(J72J) = 0.0547. 

With Hdm included, the degeneracy of the ground state 
with respect to the sign of the ordering wavevector is broken. 
Since the sign of the DM term alternates on adjacent layers, 



the direction of the wavevector alternates from layer to layer 
to give the spin structure (setting a = 0) 





sit S\ cos(Q cJ -R) 

r-l)«sin(Q cr R). 



(10) 



where now Q c / is determined by the minimum of Jq. We find 
Q c , = (0, 1/2 + t cU 0) with e cl = 0.0533. 

The classical ground state in the presence of a transverse 
magnetic field can be found easily because U(l) symmetry 
ensures that only Fourier components with Q = and Q = 
Q c ; contribute to the spin configuration. The spiral order of 
spin components within the b-c plane is preserved, and the 
spins cant towards the field direction to produce a cone state 
with 



Sr. = S I COS #0 -v. , ~ - ■ 

v (-l)"cosfl sin(Q c/ -R), 



sin 60 
) cos(Q c / • R) 



(11) 



where, measuring magnetic field in the reduced units h l = 
cf^bB 1 /S, sin 6*o = h a jh a cr . The critical field in reduced 
units is h a cr — 2 [j£ — Jq], giving B° r = 8.36 T. The same 
expression for B° r also emerges from an exact treatment of 
the quantum Hamiltonian for a single layer (see Ref. |3| and 
Sec. llV Bl ; the small difference between this value and the ex- 
perimental one 3 of 8.44 T is partly due to interlayer exchange 
(see Sec. [VJi. At higher fields, the spins are fully polarized 
along the field direction. 

Ground states in a longitudinal field are considerably more 
complex because the magnetic field breaks U(l) symmetry 
and many Fourier harmonics contribute to the spin configura- 
tion. A useful guide to the behavior one should expect is pro- 
vided by results for frustrated magnetic systems in a magnetic 
field, with single ion anisotropy rather than DM interactions. 
In that case, if anisotropy is weak, there is a first order tran- 
sition between a distorted cycloid state at low field, in which 
spins are confined to the easy plane, and an incommensurate 
cone structure with its axis along the field direction at high 
fields 

To investigate such phenomena in the problem we are con- 
cerned with, we have studied spin configurations obtained by 
minimizing the classical energy functional numerically. We 
use periodic boundary conditions with a period of over 1000 
sites in the 6-direction, and have examined many minima for 
a range of values of longitudinal fields. We find two phases 
separated by a first order transition. For h/h cr < 0.35, the 
zero-field spin spiral evolves smoothly into a distorted cycloid 
in which spins lie in the b-c plane. This state has a continuous 
degeneracy associated with phason modes. 21 For fields in the 
range 0.35 < h/h cr < 1, an incommensurate out-of-plane so- 
lution is optimal. It has a non-zero value for the chiral scalar 
K and therefore breaks Z2 symmetry. A very good approxi- 
mation to the out-of-plane numerical solution is provided by 
the expression 



(cos Qo cos(Q • R) cos 77 + sin 9q sin rj^ 
(-l) n cos6»o sin(Q • R) 
sin 9q cos 77 — cos 6q cos(Q • R) sin 77 ; 



(12) 
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In this approximation, only the Fourier components and Q 
appear, and the ordering wavevector is within a few percent of 
Q*;. Spin directions at different sites form a cone, which has 
a height S sin 8q and an axis lying in the a-c plane, tilted at 
an angle 77 to the c-direction. Moving from site to site in the 
6-direction, the spin projection onto the easy plane traces out 
an ellipse. The eccentricity of this ellipse is associated with 
a non-zero DM energy, and rj oc D for small D. A second 
ground state, related by Z2 symmetry to the first, is generated 
by the operation: Q — > Q and r\ — > —77. 

In spite of the proximity of the incommensurate wavevector 
to the commensurate value (0, 5, 0), the commensurate states 
are found to be well separated in energy from the incommen- 
surate solutions, within a classical treatment. 

The results of this classical analysis are summarized in 
Fig. Behavior in a transverse field is in qualitative agree- 
ment with the experimental findings outlined in Sec. |I] We 
delay a quantitative comparison between theory and experi- 
ment until after our discussion of the effects of quantum fluc- 
tuations in Sec. IIVI Contrastingly, observed behavior in a 
longitudinal field shows different features from the classical 
phase diagram. In particular, the state found in the field range 
2. IT < B c < 7. IT does not appear classically. 



Transverse 




FIG. 2: Phase diagram in the classical limit, with a schematic repre- 
sentation of the different phases. Transitions between the cone states 
and the ferromagnetic states are second order. The distorted cycloid 
and the tilted cone states are separated by a first order transition. 



A. Large S-expansion 

We now turn to a description of the calculations. The pro- 
cedure is standard: starting from a classical, ordered state we 
use the Holstein-Primakoff transformation to obtain a bosonic 
Hamiltoniani2i2i2^i22i22iS Considering only the quadratic 
part of this Hamiltonian, we obtain the leading quantum con- 
tribution in a l/S expansion. 



1. Transverse Field 

In a transverse magnetic field, the classical ground state is 
an incommensurately ordered spin cone with wavevector Q, 
given by Eq. \\ It . We introduce a rotating coordinate system 
in spin space, via 



s r\ 


- l 












cos(Q • R) 


-(-)"sin(Q- 






-) n sin(Q-R) 


cos(Q • R) 



R) 




(13) 



chosen so that the z axis at each site is aligned with the clas- 
sical spin direction. A central objective of this section is to 
calculate quantum corrections to classical values of the order- 
ing wavevector Q and the canting angle 9. We omit the small 
interlayer exchange J", postponing a discussion of some of 
its effects to Sec. El 

The Holstein-Primakoff transformation is 



IV. QUANTUM FLUCTUATIONS 



si = 
s R = 



2S 



2 

S-t 



(4- 



J R ) 



R 



(14) 



The classical ground states determined in Sec. 11111 provide 
a starting point for a treatment of quantum fluctuations. This 
can be approached using either a l/S expansion, or directly 
for S = I by expanding in powers of the density of reversed 
spins in a polarized background, viewing these as a dilute 
Bose gas. While the l/S expansion is uncontrolled when ap- 
plied to CS2CUCI4, it is known to produce quite accurate re- 
sults for some simpler two-dimensional S = 5 systemspA£*A&ii 
Conversely, the density of reversed spins is controlled by field 
strength and the expansion parameter is (1 — h/h cr ). It is 
worth pointing out that quantum fluctuations in CS2CUCI4 do 
not break a classical degeneracy, as is the case for the isotropic 
triangular lattice antiferromagnets in a field, 22 but are likely 
to have substantial quantitative effects on ground state proper- 
ties. 



where the boson creation and annihilation operators satisfy 



the commutation relation 
the Fourier transform 



6 R> <t>R> 



<5r,r'- Introducing 



(15) 



for a lattice of N sites, the Hamiltonian of Eq. Q takes the 
form 

H = H + H 1 + H 2 + --- , (16) 
where Ti n is proportional to S 2 ~ n l 2 and consists of products 
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of n normal-ordered boson operators. The leading terms are 
Ho = NS 2 (j£ + [j<f - j£] sin 2 6» - h a sinfl) 

Hi = y^cos0(2[J o T - jQfwO-h*) + 

H 2 = iVS ( J§ + [Jj - J§] sin 2 - ^ sin 



+ "Ml 



the orientation the plane of the cycloid. For the Hamiltonian 
we are concerned with, the DM interaction lifts the SU(2) 
symmetry and generates an excitation gap at wavevector Q, 
which becomes wider in an applied magnetic field. Recently, 
it has been shown that the spin-wave spectrum of an antiferro- 
magnet in a strong magnetic field is kinematically unstable to 
two-magnon decay. 30 Here we neglect such decay processes 
and retain only harmonic terms in the Hamiltonian. 

The ground-state energy, omitting terms O(S ) and higher, 
is then 



(17) 



where the sum on k is performed over the first Brillouin zone 
and 



1 



Ak = - [2 J k + Jq +k + JQ_ k ] - J ( 

- \ t 2 J k - J^+k - Jq-k] sin 2 
h a 



Q 



£k 



[ J o - j q] sin 2 6» + — sinfl 
- [2 J k - J^ +k - jg_ k ] cos 2 



2 [^Q+k 



(18) 



The coefficients A k and B k are even functions of k, whereas 
C k is an odd function of k. The term Hi is linear in the 
bosonic operators and vanishes if the canting angle 9 takes its 
classical value, 9q. The quadratic Hamiltonian is diagonalized 
by the Bogoliubov transformation 



0k = «k7 k + w k7l k > 
^-k = «k7 k + "k7ik> 



(19) 



where 



= l + vt = - 



^v/^B 2 



-Bv 



M k W k 



2 r>2 
k - B k 



(20) 



The diagonal form of the quadratic Hamiltonian is 

h a 



n 2 = ns{ 



[Jo - Jq] sin 2 < 



1 



k \ 7 k 7 k + j 



(21) 



where £! k = y A 2 . — B 2 . + C k is the spinwave dispersion 
relationi2i2i Setting 9 = 6q, the spectrum has a gapless mode 
at k = as a result of the U(l) symmetry. The low-lying ex- 
citations are spin oscillations within the plane of the cycloid. 
For an SU(2) symmetric Hamiltonian, a second Goldstone 
mode is present at the ordering wavevector Q of the cycloid. 
The low-lying excitations in this case involve oscillations of 



E=(H) = NS(S + 1) (J§ + [J T - Jg] sin : 
- NS(S + l/2)h a sm8 + S^2n k . 



(22) 



The ordering wavevector is to be determined by minimizing 
E with respect to Q. Following this procedure, the 1/5* cor- 
rection to E comes not only from the zero-point fluctuations 
but also from the renormalization of 9 and Q. 

Results for Q = (0, 1/2 + e, 0) to C^S" 1 ) are shown in 
Fig. |3 together with data from Ref. 3. At the critical field 
h^ r we find, in agreement with the experiment, that Q takes 
the classical value Q c ;, which is field-independent. This is a 
consequence of the fact that the ferromagnetically polarized 
state is an exact eigenstate of the Hamiltonian with vanish- 
ing zero-point energy. At lower fields, fluctuations renormal- 
ize Q, which decreases with decreasing field: the zero-field 
value of the incommensuration e = 0.21 is significantly re- 
duced from its value at the critical field. This reduction can 
be understood on the basis that zero-point energy in antifer- 
romagnets generally is lowered for states with collinear spins. 
The states we are concerned with are close to the collinear 
state with Q = (0, 1/2, 0), but have lower classical energies. 
With decreasing field, quantum fluctuations are enhanced and 
drive the incommensurate wavevector towards the commen- 
surate value. As a technical aside, we note that calculations 
are simplified by the presence of DM interactions, since with- 
out them the Goldstone mode at wavevector Q, which appears 
as h — ► 0, necessitates a self-consistent treatment of quantum 
fluctuations. 

The ordered moment is reduced from its classical value by 
quantum fluctuations. At leading order 



(S) s (Si,) 



k 



s- 



N 4^ 2 



Av 



2 Vv^T 1 ^ 



- 1 -(23) 



This is shown as a function of transverse magnetic field in 
Fig. |4] Our zero-field value of (S) = 0.25 is close to 
the result (S) = 0.266 + 0(S~ 3 ) for the isotropic trian- 
gular antiferromagnet, 17 and to results for the anisotropic 
lattice without DM interactions, obtained using the 1/5 
expansion 31 ' 32 and series expansions^ 

The canting angle 9 can be determined in two different but 
equivalent ways. Classically, the condition 9 = 9 ensures 
both that Hi = and that (Ho) is at a minimum. The leading 
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FIG. 3: The incommensuration e as a function of transverse mag- 
netic field strength B a . Full line: result from 1/5* expansion. Long 
dashed line : result from the classical theory. Short dashed line: lin- 
ear variation of e with B a , from calculation for dilute Bose gas of 
spin flips. Symbols are the experimental results taken from Fig 3c 
of Ref. 3 taken at T = 0.20K: A from magnetic Bragg peaks at 
Q = (0, 1.5 - e, 0), and □ from peaks at Q = (0, 0.5 - e, 1). 



1/S correction can be determined similarly. First, normal or- 
dering of H3, expressed in terms of 7^ and 7 k , yields a term 
linear in boson operators, which should be combined with Hi : 
the combination vanishes when 9 takes its ground state value. 
Second, and alternatively, one can minimize (Ho + H2) with 
respect to 9. In this way we find 



sin c/q 



(24) 




where A k and B k should be evaluated at 6q. Because quantum 
fluctuations are suppressed as the critical field is approached, 
9 — ► 6*o as h a — > h" r . As seen in the inset of Fig. |4j the 
quantum corrections to sin 9 are small and Eq. [24] is nearly 
equal to the unrenormalized function sin 9q(~ h a jh a cr ). 

Combining results for the ordered moment and the canting 
angle, the magnetization is given by 



m 



9 

N 



52{S&=g a iiB{Stism6, (25) 



R 



which yields 



1 



B k A k -B k 



2SN ^ B 

k 



A k + B k 



(26) 

for h a < h a cr , and m a = g a fi B S for h a > h a cr . The 1/5 cor- 
rection on the right-hand side of Eq. J26b has a dependence on 
the magnetic field through the values of A k and B k (which 
again should be evaluated at 9 ). The departure of the mag- 
netization curve from the simple linear dependence expected 



FIG. 4: Ordered moment as a function of transverse magnetic field 
strength. Dashed line: classical behavior; full line: with leading 
quantum corrections, from Eq. I23i . In the inset, Ea. l24l is plotted as 
a function of the transverse field in full line, whereas sin #0 is plotted 
is dashed line. 



classically is hence a consequence of zero-point fluctuations. 
To understand the sign of this departure, it is useful to recall 
that the ground-state energy E(B a ) as a function of field is 
related to the magnetization m a (B a ) via 



E(0) - E{B a cr ) = N / m a (B a )dB a 



(27) 



Now, since fluctuations reduce E(Q) below its classical value 
but do not contribute to E(B® r ), the fluctuation contribution 
to Eq. d27i is negative. Supposing this to be true not only for 
the integral but also for the integrand at all B a , it is natural 
to expect the magnetization curve at finite S to lie below the 
classical one for all B a < B^ r . A comparison between our 
results and experimental data, presented in Fig. |5] shows a 
very good agreement. 




FIG. 5: Magnetization as a function of transverse magnetic field 
strength. The dashed line shows the result of classical theory. The 
thick, full line includes the 1/5 correction from Eq. 1261 . The thin, 
full line gives experimental data, 34 measured at T = 60mk. Inset: 
calculated ground state energy as a function of transverse magnetic 
field. The dashed and thick lines are the results of classical theory 
and the 1/5 correction, Eq. 1221 . respectively. 
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The component S T of the ordered moment in the plane per- 
pendicular to the applied field can also be evaluated within the 
1 / S expansion. Defining it by 



S T = |<SR-a(S£) 



(28) 



we find S T = (S) cos 9. Using Eqs. (I23> and d24> . we obtain 
at first order in 1/5 



S 1 = Scosflo 



sec 2 9 
2SN 



E 



-4, 



- 1 



tan 2 O sr-^ Bk 
2.^ ~rT 



2SN 



-Bq V ^4k + -Bk 



(29) 



2. Longitudinal Field 



In longitudinal field, calculations of fluctuation effects us- 
ing the 1/5 expansion are complicated by the fact that the 
classical ground state contains many Fourier harmonics. At 
low transverse fields, the classical ground state consists of 
a distorted cycloid in which spins lie within the b-c plane, 
as described in Sec. |lll| In this regime we write Sr = 
5(0,cos<^R,sin(/iR,) and consider the leading anharmonic 
distortion to the cycloid structurep2i2& 



f> R = Q R + /3cosQ R, 



(30) 



Results are presented in Fig. |6] While classical theory 
gives S T oc (1 — (B a /B^ r ) 2 ) 1 ^ 2 , fluctuations generate a non- 
monotonic dependence of S T on B a at low fields. This be- 
havior can be understood on the basis that polarization of the 
spins with increasing applied field has the effect of reducing 
the phase space available for quantum fluctuations, and hence 
increases order. Experimental data are also shown in Fig. [6] 
since the absolute scale for S T has not been determined, we 
scale the data to fit theory at high fields. The result of the 
1/5 expansion compares favorably to the experimental data, 
which also shows that at low field the perpendicular ordered 
moment increases with increasing field . 



0.5 ■ 



. 4 



where for concreteness, the field is taken to be along the c 
axis. The distortion of the cycloid is parameterized by j3: its 
value, determined by minimizing the classical energy, is 



J 2Q 



J T - 2 J T 



(31) 



We consider quantum fluctuations about this classical state, 
using the Holstein-Primakoff transformation and omitting 
terms 0([h c ] 3 ) and 0(5°) to obtain the Hamiltonian 



H 
c/3 



0.3 



. 2 : 



a n 



(32) 



o.i 



with 



FIG. 6: Component of ordered moment S in the plane perpendic- 
ular to the field direction, as a function of transverse magnetic field 
strength. Dashed line: classical theory. Full line: result including 
1/S corrections, from Eq. I29i . Symbols (□): experimental data at 
T < 0.1K (Taken from Fig. 3b of Ref.0) 



and 



Ho 



NS 2 



(h 



CY2 



J 2Q 



JT 



2^ 



(33) 



H 2 = NSJ% + SJ2 Uk(Vl> ; 



iSd \ 



k+Q 



-k^k-Q 



h.c 



[£k- 



-k^k 



J 2Q 



J, 



k-Q 



-k^k-Q 



h.c. 



(34) 
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Here 



A k = J [2 Jk + J^ +k + j£_ k ] -J% + (3 



B L - J [ 2 ^k - JQ+k - ^Q-k] - & 



2 | J 2Q+k ^ J 2Q-k 



16 



2 I J 2Q+k T J 2Q-k ' ,7 k ' "'-k 
16 



tT , jT 
J Q+k + J Q-k 



rT , tT < 
J Q+k + J Q-k 



( ^3Q-k + ^ T 3Q+k ^ J T 2 Q+k ^Q+k ^-Q-k + ^J-k + 2 </lq +k ) 



(35) 

(36) 
(37) 



The higher Fourier harmonics in the classical ground state 
scatter spin fluctuations with a momentum transfer which is a 
multiple of Q. The presence of these scattering terms, propor- 
tional to (3, in the quadratic spinwave Hamiltonian means that 
the dispersion relation is determined by an infinite set of cou- 
pled equations. 35 Since our calculation is anyway restricted to 
small h c oc (3, we treat these coupled equations to C(/3 2 ) in 
a calculation of the ground state energy. More specifically, 
it is convenient first to perform a Bogoliubov transformation 
to diagonalize the momentum conserving terms, and then to 
evaluate the contribution from the terms scattering by ±Q us- 
ing perturbation theory. We find 



E = NS(S- 



NS 2 (h c f 



["^2Q + Jq 



^k + ^Q-k 



With Cl' k = 



(BTp and 



(38) 
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FIG. 7: The incommensuration e as a function of longitudinal field 
strength. Solid line: result from 1/S expansion, from Eq. <38> . 
Dashed line: result from the classical result. Symbols are the experi- 
mental results taken from Ref.|2, Fig. le: A from magnetic peaks at 
Q = (2, 1/2 + e, 0) and □ from peaks at Q = (1, 1/2 + e, 0). Ex- 
perimentally, a distorted cycloid occurs for B c < 1.4T and a second 
incommensurate phase occurs for fields in the range 1.4T < B c < 
2. IT. 



4 = - g [ L k« 



Q-k u Q~k) 



^q)««q-k + u k w Q-k)L 



(39) 



where it k and v'^ are given by Eq. J20i after substituting A^ 
and for A k and f? k . 

The ordering wavevector and its dependence on field can 
be calculated by minimizing Eq. (I38> with respect to Q. It is 
interesting to note that, in contrast to the case for a transverse 
field, the ordering wavevector in a longitudinal field is depen- 
dent on field strength even at the classical level. It increases 
with field and this trend is reinforced by the quantum fluctua- 
tions. Results are shown in Fig. [7] together with experimental 
data. The observed increase of Q with field is much faster 
than the calculated one; the origin of this discrepancy is not 
understood. 

Now we turn to the case of stronger longitudinal fields. In 
the field range 0.35 < h c /h c cr < 1, the classical treatment 
described in Sec. |in| gives the tilted cone of Eq. (I12> as the 
ground state. The tilting angle r] of the cone axis from the field 
direction is given approximately by tan 77 w (D/h c cr )(l — 
(/z° r //i c ) 2 ) and is less than one degree for h c > 0.9h^. r . Ex- 
perimentally, an incommensurately ordered state has recently 
been observed 5 for 7. IT < B c < B c cr . With this in mind, we 
approximate the classical ground state in this field range by 



setting 77 = in Eq. (11 2t and use the 1/S expansion to study 
the effects of quantum fluctuations. Following a procedure 
similar to the one described for a transverse field, we have cal- 
culated the ordering wavevector as a function of field. While 
the observed phase has not so far been fully characterized, its 
ordering wavevector has been measured as a function of field 
strength. From our classical calculation, we expect the state 
to be a tilted cone. Our results for the ordering wavevector are 
compared with experimental data in Fig. [8] Calculated and 
observed values of the ordering wavevector vary in the same 
way with field, but there is an offset between the two which 
remains a puzzle. In the following subsection dlVB> . we ob- 
tain results which are essentially exact close to B c cr . Since the 
discrepancy remains, we conclude that the value of the order- 
ing wavevector is influenced by interactions not included in 
the model Hamiltonian of Eq. Q. 



B. Dilute Bose Gas 

An alternative to the 1/S expansion can be motivated by 
noting that fully polarized states are exact eigenstates of the 
Hamiltonian. The absence of quantum fluctuations suggests 
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FIG. 8: The incommensuration e as a function of longitudinal mag- 
netic field strength B c . Solid line: result from 1/5 expansion. 
Dashed line: result from calculation for dilute Bose gas of reversed 
spins. Symbols (□): experimental data 5 (no incommensurate order- 
ing is observed for 2.1 T< B c <7.1 T). 



a systematic expansion in powers of the density of reversed 
spins, 36 37 equivalent to an expansion in powers of (1— hj h cr ) . 
In this approach reversed spins constitute a dilute gas of 
bosons with hard-core repulsion. 



in a transverse field. Magnetic order of in-plane spin com- 
ponents below the critical field translates, using Eq. j40> . to 
formation of a Bose condensate for [i > 0. We introduce the 
order parameter tpq and shift the boson annihilation operator 
by a constant 



4>Q 



NlpQ + (j)Q , 



(45) 



where the ordering wavevector Q need not take the classical 
value Q c ;. Minimization of the ground-state energy is equiva- 
lent to requiring (4>q) = 0. Working in the low density limit, 
the scattering amplitude between bosons is given by an effec- 
tive interaction potential r q (k, k') which results from sum- 
ming ladder diagrams. It satisfies the integral equation 



r q (k,k')=y * ^ 



Vr 



q-q' 



q' k+q' 



t k'-q' 



— iV(k,k'). 

e k' 

(46) 

The ground state energy, including the leading interaction 
effects at low density, is 



E {2) =N 



+ (3 m) I^q l 2 + ^r (Q, Q) \iPq\ 4 



5 £(* 



*i0 



(47) 



1. Transverse Field 



where £ k = ^/Ff 



~GZ + N k , and 



We introduce boson creation and annihilation operators, 
and <f) R to represent spins, setting, for a transverse field, 

2 _ ^R^R' 



4 



s R = 



R — 



4 



(40) 
takes 



with the constraint that the particle number riR = </>r_0 r 
only the values and 1 . This is imposed by introducing an on- 
site interaction U and taking the limit U — ► oo. 

The Hamiltonian for a single layer (with, for definiteness 
the layer index n chosen to be even) is 

k 



n = n 



i 

2N 



k,k',q 



VJk,k') 



where 



c k 



J k 

h a 



- h a 



(41) 

(42) 

(43) 
(44) 

Standard techniques 38 developed for the interacting Bose 
gas can be applied to treat this Hamiltonian. For h a > h a cr 
the spin system is fully polarized. Equivalently, for [i < 
the ground state is the boson vacuum. Conservation of boson 
number follows from U(l) symmetry of the spin Hamiltonian 



Vq(k,k') = 2J q + 2{7. 



Q+k 



fc Q-k 



r k (Q,Q + k) 



T_ k (Q + k, Q) + r (Q, Q + k) + r (Q + k, Q) 



Gk = |V>Ql 2 r k (Q,Q) 



N k = 



(48) 



The condition (<^>q) = yields an expression for the order 
parameter 



r (Q,Q) 

Substituting this into Eq. ( 1471 . the ground-state energy is 



(49) 



£< 2 



N 



4 



h a 

T ~ 2r (Q,Q) 



2 ^ 

k 



(-E k - F k ) 



As in the 1/5 calculation, the ordering wavevector can be 
determined as a function of field by minimizing with re- 
spect to Q. Our focus here, however is on exact results close 
to the critical field. At h a — h a cr we find Q = Q c ;. In addi- 
tion, we obtain 



dQ 

dh a 



h a =h 



"J 4 



d 3 E^ 



dQ 2 k dfi ) dQkdii* 



and hence 
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6>Q 



h a =hi 



4 



dQ k 2 



dQ, 



■In 



r (Q,Q) N 



1 E- 

N £—1 



Q+k 



c Q-k 



r k (Q,Q) 
r (Q,Q) 



(50) 



Q=Qc 



A potential difficulty arises at this point because interac- 
tions are marginally irrelevant at the critical point of the two- 
dimensional Bose gas: 39 in consequence, ro(Q c /,Q c ;) van- 
ishes for an isolated layer. It is therefore essential to in- 
clude interlayer exchange J" in the calculation of the ver- 
tex. Evaluating Eq. (15 0> numerically, we find h^ r dQ/dh a = 
(0,0.0911,0). This result is displayed in Fig. [3] it is sim- 
ilar to that given by the l/S expansion, indicating that the 
linear spinwave theory captures the effects of quantum fluc- 
tuations quite accurately in this system. Both approaches are 
in good agreement with experiment, especially close to the 
critical field. 



2. Longitudinal Field 

A similar procedure can be followed for the system in longi- 
tudinal field (chosen along the c axis, without loss of general- 
ity for the Hamiltonian of Eq. Q). With this field orientation, 
the expressions for spin operators in terms of Bose operators 



R 



^R 



R 



S 



R 



J R' 

4 



(51) 



The Hamiltonian for a single layer (again taking the layer in- 
dex n to be even) is 



n 



j - h c 



k 

2 7^ E ^ k + D ^ + h - c ) 

k,k' 

l - V q (k,kO<4 +q 4,_ t A,<£ k) (52) 



-qVk'Vk' 



k,k',q 



where 



4 - H = Jk- Jo + h c /2 



(53) 



The presence of a term cubic in boson operators consider- 
ably complicates the analysis, since with it, particle number is 
not conserved. Its appearance reflects the fact that a longitu- 
dinal field breaks U(l) symmetry as discussed in Sec. [n| In 
general, the particle number (or longitudinal magnetization) is 
not conserved (except for the boson vacuum - the ferromag- 
netic state - which is an exact eigenstate of the Hamiltonian). 
The remaining, Z2 symmetry is invariant under the canonical 
transformation </> k — > — 0_ k . 



We note in passing that the cubic term does not result in 
a first order transition from the fully polarized state as /i is 
varied, because momentum conservation precludes contribu- 
tions involving only the ordering fields <pQ and 4>-q in a 
Landau-Ginzburg description. An ordered state is therefore 
brought about by the closing of the single-particle excitation 
gap, yielding a second order phase transition. The universal- 
ity class associated with this quantum phase transition must 
take into account the extra Z2 symmetry of the Hamiltonian. 
The low energy action is described by a Z2 © U(l) symme- 
try model. This multicritical transition found in longitudinal 
field is to be contrasted with the ordinary XY quantum phase 
transition found in transverse fields 

While for a transverse field the ordering wavevector can be 
found simply from the quadratic part of the boson Hamilto- 
nian, this is not so for a longitudinal field. In that case, because 
particle number is not conserved, the quasiparticle spectrum 
is renormalized by quantum fluctuations, even at the critical 
point. It is interesting to check whether a renormalization of 
this kind may be responsible for the discrepancy between the- 
ory and experiment shown in Fig. [8] The critical field and the 
ordering wavevector are determined from the values of h c and 
Q for which the one-particle Green function has a pole at zero 
energy, by solving 



G(Q,£ Q = 0) =0. 



(54) 



In absence of DM interactions, the one-particle Green func- 
tion at and above the critical field is given exactly at zero 
temperature by the expression for a non-interacting system, 
G° (k, iuj) = (iu> — e£ + Since D -C J, we evaluate 

the leading contribution to the self-energy, 



E(k,i W ) = -L^ r i 

2iv ■* — ' tuj - 
q 



Di 



2[i 



O (D 4 ) 



(55) 



yielding a renormalized quasiparticle spectrum ~ e£ — //-+ 
E(k, e k — /i) in the symmetric phase. The ordering wavevectoi 
can be found by solving Dyson's equation at the critical field, 



G° (Q,0)- 1 -S(Q,0) = 0, 



(56) 



which gives Q = Q* ; + (0, 0.00025, 0). This minute quantum 
correction at the critical field is nearly two orders of magni- 
tude too small to explain the discrepancy between Q*, and 
the experimental ordering wavevector illustrated in Fig. [8] We 
conclude that there are further anisotropic interactions present 
in the system but not captured by the Hamiltonian of Eq. Q. 
Additional evidence for this is provided by the fact that the 
experimental phase diagram in a longitudinal field depends 
on field orientation within the b-c plane. 5 
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It is interesting to note that at the critical point in a lon- 
gitudinal field, in contrast to behavior for a transverse field, 
order is possible at two wavevectors, ±Q. Cone states break 
spontaneously the Ising symmetry, with condensation either 
at Q or at — Q. An alternative possibility is the simultaneous 
condensation of magnons at both wavevectors, forming a fan 
phase. Competition between the fan and cone phases is deter- 
mined by the interaction between magnons. A straightforward 
calculation shows that the cone phase is favored if 

r„(Q,Q) <r (Q,-Q) + r 2Q (Q,-Q), (57) 

while the fan phase is preferred if the inequality is reversed. 
Evaluating the vertices numerically, we find that, although 
quantum fluctuations renormalize interactions they do not 
modify the character of the ground state found from the clas- 
sical calculation, and the cone state is favored. 



V. INTERLAYER COUPLING 

To this point, we have omitted the interlayer coupling J" 
(except where it was essential, in order to obtain a non-zero 
value for the interaction vertex r q (k, k')). It is relatively 
weak (</"/ J ~ 0.05), though crucial in stabilizing long-range 
order against thermal fluctuations. As well as being small, it 
is also frustrated by DM interactions, because the sign of the 
DM interactions alternates between layers (see Eq. l|3}). The 
frustration introduces distortions in the cone states, which we 
discuss in this section. 

More specifically, considering zero field for simplicity, 
the classical ground state in the absence of interlayer cou- 
pling consists of a spin spiral with wavevector +Q in layers 
with even index n, and wavevector Q in odd layers, as in 
Eq. ( I10i . By contrast, for a system with antiferromagnetic 
interlayer exchange but no DM interactions, the ground state 
consists of spin spirals with the same wavevector (say +Q) 
in every layer, and with alternating phases a in even and odd 
layers, so that 





■>)( = | cos(Q ■ R + mi) 
sin(Q • R + wr) 



(58) 



With both interlayer exchange and DM interactions, their 
competition results in a ground state which is a superposition 
of the two structures. 3 - In the presence of a transverse field, 
spins lie on an elliptical cone around the field direction with 



Sr = (- 



(-1) 



l) n Si cosQ 
Si sin Q 

■? (R) 
l)"S 6 cosQ ■ 



<S a (R) 
R+(-l) n S 2 cos(-Q- 
R + S 2 (sin-Q-R) 

R 



R) 



(59) 



S c sin Q • R 



where U(l) symmetry has been broken by selecting the b- 
component of the spin to alternate on adjacent layers. The 
eccentricity I = (S% + S 2 )/(Si - S 2 ) = S b /S c of the cone 



is a measure of the mixing of the two spin spirals at wavevec- 
tors ±Q. 

Experimentally, this ratio can be measured by determining 
the relative intensity of two magnetic Bragg peaks associated 
with the spin ordering, 5 and we focus on its field dependence. 
The mixing between the two spin structures is observed to be 
particularly strong near the critical field, where it can be cal- 
culated using linear spinwave theory. 3 More generally, we find 
the field-dependence of / by minimizing the classical energy 
over states which are parameterized as in Eq. J59I . Results 
(obtained numerically) compare well with experimental data, 
as shown in Fig. [9] Mixing is small (/ ~ 1) in zero field, but 
rises rapidly near the critical field to reach the value I = 1.52 
atB a = B a cr . 

To gain insight into these results, it is useful to consider 
behavior close to the critical field, and expand in powers of 
the small in-plane spin components, Sb and S c . Following 
this procedure we obtain to quartic order the energy 



E 



N 



S 2 (Jo + J" - h a ) 



+ ( Jq - Jo - J" + h a /2) (S 2 + S 2 ) /2 
+ D Q S b S c + J" (S 2 S 2 ) /2 + (S 2 - S 2 C ) 



h a 
MS 2 



(60) 



This should be minimized with respect to Sb and S c . It is 
convenient to change variables, writing S b = r sin \ and S c = 
r cos %, so that 



N 



S 2 (Jo + J" - h a ) 



( Jq - Jo - J" + h a /2 + D Q sin 2 X + J" cos 2 X ) 



4 / J2Q - Jo 2 o 

r I — _ — cos 2x 



32S 2 



6AS 2 



1 + 2 cos 4 x + 2 sin 4 x 



(61) 



The eccentricity is then I = tanx. As the crit- 
ical field is approached from below, r — > + and 
X is determined solely by the quadratic term, yielding 
tanx = Dq/(J" - ^J{J") 2 + (D Q ) 2 ) = 1.52, as reported 
previously. 3 Note that the interlayer exchange modifies the 
previous estimate for the critical field (see Sec. IIIII to h a cr — 

2[Jo - Jq + J" + V(J") 2 + {Dq)' 2 }, giving B a cr = 8.51T. 
With reducing field, r increases and / is determined partly by 
the x dependence of the quartic term, which is minimum at 
X = 7r/4 + mir/2 (taking h a ^> J 2 q — Jo). The quartic term 
hence favors |/| = 1 and dominates as h a is reduced below 
h a 



VI. CONCLUSIONS 

In summary, we have presented a detailed investigation 
of ground-state properties of an anisotropic triangular lattice 
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FIG. 9: Eccentricity of ellipse / = Sb/S c as a function of the trans- 
verse magnetic field strength. Line: result of classical analysis. Sym- 
bols (□): experimental data!" 



Focussing on the ordering wavevector of the incommensurate 
phases, spinwave theory gives only a mediocre account of its 
behavior at low field, while neither spinwave theory nor cal- 
culations for the almost polarized system can explain its value 
close to saturation. In this connection, it is worth emphasizing 
that the model Hamiltonian we have used must in fact omit 
some residual interactions which are of importance, since it 
is invariant under rotations of the magnetic field about the a- 
direction, while the observed phase diagram does not have ex- 
actly this symmetry^ 

For the future, the nature of the ground state in a longitudi- 
nal field at intermediate field strengths remains an intriguing 
problem, which we intend to address elsewhere^ 

The work was supported by EPSRC under Grant 
GR/R83712/01 (JTC) and Grant GR/R76714/01 (RC). 



APPENDIX A: LADDER DIAGRAM SUMMATION 



antiferromagnet with Dzyaloshinskii-Moriya interactions, fo- 
cussing on behavior in an applied magnetic field and its depen- 
dence on field direction. We have supplemented calculations 
for a classical model with two approaches to quantum fluctua- 
tions: one using linear spinwave theory, and the other treating 
reversed spins in an almost polarized state as a dilute Bose gas. 
We have compared our calculations with experimental data on 
Cs2CuCl^i*2iii^i The outcome of this comparison depends 
strikingly on field direction. For a transverse field, theory is in 
qualitative, and on many points quantitative agreement with 
experiment. For a longitudinal field, central aspects of the 
low-temperature phase diagram remain to be understood. 

In more detail, for a transverse field the classical model 
yields the observed incommensurate cone state with a field- 
dependent canting angle. Interlayer exchange interactions in- 
fluence this ordering in ways that are also well-described by 
classical theory. However, to account for the measured field 
dependence of ordering wavevector, the magnetization and the 
local ordered moment, it is necessary to include the effect of 
quantum fluctuations. Linear spinwave theory gives quite ac- 
curate results for the magnetization and qualitatively correct 
behavior for the other two quantities; the large reduction at 
low fields of the ordered moment below its classical value 
demonstrates the importance of fluctuations. In addition, cal- 
culations for the almost polarized system fit observations of 
the ordering wavevector very well, as they should since that 
aspect of the theory is asymptotically exact. 

By contrast, for a longitudinal field there are clear differ- 
ences between the phase diagram of the classical model and 
experiment. In this case, classical theory yields a distorted 
cycloid as the ground state at low fields, separated by a first- 
order phase transition from a tilted cone state at higher fields. 
Experimentally, most of the high-field region is occupied by 
a third phase, in which no magnetic Bragg peaks have been 
reported, 3 although incommensurate order has recently been 
observed in a narrow field window below the saturation field. 5 



In this appendix we show how we solve numerically 
the integral equation for the effective interaction potential, 
r q (k, k'). We recall Eq. g6) 



(Al) 



where for clarity we have omitted the variables k, k' and in- 
troduced 



1 



Ek4 



Ck'-c 



(A2) 



The bare interaction Vq arises as the Fourier transform of an 
interaction in real space, in the form Vq = J2r ^R exp(— iq- 
R). Crucially, since the interaction is short range, only a small 
set of coefficients ^4r are non-zero. In turn, this implies that 
r q can also be expressed using a small number of Fourier co- 
efficients, as follows. Define the parameters £?r through the 
equation T q = J2r Br exp(— iq • R). Then from Eq. iAli 
we obtain 



Br = A 



R 



1-]TM R , R ,5 R , 

R' 



(A3) 



with M R r/ 



X) q n q e lq ^ R A simple consequence 
of Eq. JA3t is that if /1r = for a given R, then _Br = as 
well. From Eq. iA3\ we find 



B R = 



R' 



M 



R',R 



(A4) 



Since the Hamiltonian has only nearest-neighbor interactions 
on a stacked triangular lattice, the matrix we must invert has 
9x9 elements. These can be evaluated numerically. 
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